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Dissipation coefficients are calculated in the adiabatic, near thermal equilibrium regime for a 
large class of renormalizable interaction configurations involving a two-stage mechanism, where a 
background scalar field is coupled to heavy intermediate scalar or fermion fields which in turn are 
coupled to light scalar or fermion radiation fields. These interactions are typical of warm inflation 
microscopic model building. Two perturbative regimes are shown where well deflned approxima- 
tions for the spectral functions apply. One regime is at high temperature, when the masses of both 
intermediate and radiation fields are less than the temperature scale and where the poles of the 
spectral functions dominate. The other regime is at low temperature, when the intermediate field 
, masses are much bigger than the temperature and where the low energy and low three-momentum 

' regime dominate the spectral functions. The dissipation coefficients in these two regimes are derived. 

' However, due to resummation issues for the high temperature case, only phenomenological approx- 

Qj^, imate estimates are provided for the dissipation in this regime. In the low temperature case, higher 
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loop contributions are suppressed and so no resummation is necessary. In addition to infiationary 
cosmology, the application of our results to cosmological phase transitions is also discussed. 
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I. INTRODUCTION 



Many physical processes involving adiabatic dissipative effects near thermal equilibrium occur in relativistic systems 
of bosonic and fermionic quantum fields. In the evolution of the early universe, several phase transition are believed 
, to have occurred, with Grand Unified, Electroweak, quantum cromodynamics (QCD), and chiral amongst the most 
^-H ■ commonly studied, and for both first and second order type transitions. Many physical implications are associated 
with these transitions including baryogenesis, defect formation, magnetic field formation, and relic particle production. 
In all such cases, interactions of the order parameter, typically a scalar field, with other fields will lead to various 
^) ' dissipative effects, and these effects, though not always accounted for, can have significant influence on physical 
processes. 

An example of a phase transition in which dissipation has significant influence and even changes the basic dynamical 
nature of the process has been seen in the inflationary universe scenario. This scenario involves the evolution of a 
scalar field which during the process of changing phases of the system also induces inflationary expansion of the 
H universe. In the standard picture of inflation [Ij, the inflaton field is modelled such that the interactions with other 
fields play no role during inflationary expansion. This leads to a thermodynamic supercooled phase in the universe 
during inflation. After this, to end inflation and put the universe back into a radiation dominated phase, the inflaton 
enters a period of coherent oscillations, where interactions with other fields now become important to reheat the 
universe with particles An alternative to this cold inflation picture is to consider dynamics where the inflaton 
interacts with other fields all during inflation. In such a case, dissipative effects become important and the dynamics 
of inflation is completely altered. In this warm inflationary dynamics , radiation production occurs concurrent with 
inflationary expansion, and the seeds of density perturbation are thermal in origin 4-6], thus signiflcantly altering 
the observable signatures of inflation compared to cold inflation. 

Our primary motivation for the dissipative coefficients calculated in this paper is for their application to warm 
inflationary cosmology. However, our results may also have a broader range of applications in different contexts 
involving the nonequilibrium dynamics of flelds, and we will briefly discuss these possibilities later in the paper. 
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Typical warm inflation models involve a two-stage field interaction mechanism, where a scalar field is coupled to 
another intermediate field which in turn is coupled to yet another decay field. The intermediate field (which will be 
called the "catalyst" field) is a heavy field (in this paper we consider it either a scalar or a fermion) that can decay 
into light fields (which will be called "radiation" fields and will here again be either scalars and/or fermions). These 
types of interactions, although naturally emerge in the hierarchy of interactions, were first realized for calculations of 
dissipation coefficients from model building applications to warm infiation and the interaction structure was generically 
termed as the two-stage mechanism 7] (for recent reviews, please see @, Q). These two-stage interactions, although 
more complicated than the direct interactions, where the relevant field modes would directly decay in light radiation 
particles, offer some simplifying features. For one thing, in the regime where the intermediate field mass is heavy 
relative to the temperature whereas the radiation fields are light, what we term the low-temperature regime, we 
can show that leading order contributions suffice in general to give a reliable estimate for the dissipation coefficient. 
This may not be the case in the regime where both the intermediate and final field masses are small relative to the 
temperature, what we term the high-temperature regime, where a different coupling governs the vertex compared to 
the decay width and the appearance of infrared and near on-shell divergences preclude the need for resummation of 
higher order terms, as we will detail later on. 

In this paper we compute dissipation coefficients coming from intermediate fields, either scalars or fermions, in 
two-stage field interacting models, in the adiabatic, near thermal equilibrium regime. Within this work, we develop a 
very new approximation which we call the low-momentum approximation. We show in this paper, both analytically 
and numerically that this low-momentum approximation is valid for the low-temperature regime, where a heavy mass 
field is involved in mitigating the interaction between the background field and the light fields. This low-momentum 
approximation was first applied in . In this paper we will further extend those earlier results and give a detailed 

analysis of the possible different temperature regimes that can be defined in terms of the different energy scales 
involved. 

To contrast the analysis in this paper of this new low-momentum approximation, there is another very common ap- 
proximation used in the literature, which picks out the poles of the propagators as dominating the integral expressions 
of dissipative coefficients. We will refer to this as the pole approximation. One key result we will show in this paper 
is for these sorts of two-stage processes where the scalar background field interacts with the light radiation fields via 
heavy catalyst fields, in the low-temperature regime, defined as where the heavy catalyst field masses are much bigger 
than the temperature scale, this new low-momentum approximation dominates over the pole approximation. In par- 
ticular, if one extended the pole approximation down into this low-temperature regime, one would find exponentially 
damped behavior for the dissipative coefficients. However, since the low-momentum approximation dominates this 
effect, the result is that the dissipative coefficients damp much slower, as a power law with respect to temperature, 
and so are much bigger in the regime of low-temperature, as compared to the results that would be obtained from 
the pole approximation. 

This paper focuses only on renormalizable two-stage interactions between scalar and fermion fields. A thorough 
analysis of dissipation coefficients for many of the applications stated above, also require treating gauge fields, which is 
beyond the scope of this paper, but we plan to examine them in a future paper. Nevertheless, even in gauge theories, 
there are interactions of the sorts computed in this paper, which up to now have not been treated in the literature. 

The paper is organized as follows. In Sect. |lT]we briefiy review the physics and definitions leading to the appearance 
of dissipative terms in the effective equations of motion for background fields. In Sect. Illll we specify all the interaction 
Lagrangian densities which we will treat in the paper. Particular attention is given to the interactions relevant to 
the two-stage mechanism Q , where dissipative effects are determined by the decay of a heavy intermediate field into 
light fields, with the former coupled to a relevant background field, whose dynamics we are interested in. The formal 
expressions for calculation of the dissipative coefficients, within the model interactions and decay forms we considered 
in our derivations, are given in Sect. IIVI In Sect. |V]we give a physical interpretation of the dissipation coefficients 
by showing their connection with particle production. In Sect. IVII the spectral functions, which are the fundamental 
quantities needed for calculating the dissipation coefficients, are computed. The results for dissipation coefficients in 
the very low, low and high-temperature regimes are presented in Sect. IVIII Various applications of these coefficients 
to inflation and phase transitions in general are examined in Sect. IVIIII Section HXl gives our summary and concluding 
remarks. Two Appendices are also included to give technical details. In Appendix[X]we show possible model building 
realizations of the type of interactions we consider in this paper. In Appendix |B] we give some of the details for 
the calculations of the imaginary parts of the relevant self-energy contributions considered in this paper and the 
corresponding particles decay widths. 
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II. DISSIPATION COEFFICIENTS 



In this section we describe how dissipation coefficients are defined in quantum field theory. We start by discussing 
the definition and derivation of the dissipation coefficient in a generic interacting quantum field theory model. 

A system, displaced from a state of equilibrium, when interacting with an environment, is expected to exhibit non- 
unitary evolution. The system dynamics will typically experience dissipative effects as it returns to its equilibrium 
state. A common example involving quantum fields is the case of some background scalar field (/? that is interacting 
with other fields, with Lp having an amplitude that is initially displaced from equilibrium. Let us describe this field 
theory model as having standard kinetic terms and with a potential given by 



V{^,X)=V{^)+Vn,t{^,X) , (2.1) 

where X represent any other field or degree of freedom that is coupled to the background field and VintifiX) gives 
how this background field is coupled to this other field. A proper study of the evolution of this background field can 
be performed in the context of the in-in, or the closed-time path (CTP) functional formalism (for an introduction, 
see for example [iT]). By integrating over the X field, a nonlocal effective equation of motion for ip can be derived. 
If we write the interaction term in Eq. ()2.ip as 



yi„t(^,X) = /(^)5(X), (2.2) 

then an ensemble averaged effective equation of motion for the system, described by the background field has the 
generic form up to first order in the nontrivial nonlocal terms (see also Ref. 12] for a recent review) 

^'^^1^ - J " ''^ [^(^("'^^ ^ •^(^("^^1 = " ' ('-^^ 

where £(,ff,r['^] is the renormalized effective Lagrangian density for (p and (• • •) are ensemble averages with respect to 
an equilibrium (quantum or thermal) state. These ensemble averages over the fields can also be conveniently expressed 
in terms of the causal two-point Green's functions for the fields X, as shown e.g. in Refs. (tI. [sI. [l3l - [T6| . with the form 
of the nonlocal term in Eq. (|2.3p dependent on the form for the interaction term Vint{^, X). A simple case is for 
example when we have a bi-quadratic interaction of the scalar field if, assumed homogeneous, ip = ipit), with some 
other scalar field represented by X, 

V,ntiv,X)^^^^X^ . (2.4) 
In this case, the nonlocal term in Eq. (|2.3I) becomes [l^ 



/ '^''^'^(^ " ifM^')) - /(^(^))] {[g{X{x)),g{X{x'))]) 
^2gt^[t) f Jt'[p{t'f~vm J |^Im[G++(p,t,i')]L' ' (2-5) 

where G~\^'^ {p,t,t') is the spatial Fourier transformed causal two-point propagator for the X field. 

The non-local term in the effective equation of motion represents a transfer of energy from the ip field into radiation, 
i.e. it is a dissipative effect. Such terms can be localized when there is a separation of timescales in the system (for 
discussion about the validity of the localization approximation for the nonlocal term and region of parameters where 
this can be achieved in specific models, see e.g. Refs. [I,[l3l)- Suppose, for example, that the self-energy introduces 
a response timescale r. If tp is slowly varying on the response timescale t, 

^«r-i, (2.6) 

which is typically referred to as the adiabatic approximation, then we can use a simple Taylor expansion and write 

fMt')) fMt)) = it' - 1) + ■■■ (2.7) 



4 



The effective equation of motion for Lp (assumed homogeneous) including the linear dissipative terms then becomes 

^ + T^+^%M^O, (2.8) 

dip 

where VeS^vi'^) is the renormalized effective potential and T is the dissipation coefficient defined as 

T = j (fx'Y.R{x,x'){t' -t), (2.9) 



where S_R(a;, x') is a retarded correlation that depends on the specific form of the interaction Vint(</?, X) and is defined 
by 



T;b{x,x') = -i 



dm 

dip 



e{t-t'){[g{X{x)),g{X{x'))\) . (2.10) 



III. INTERACTION TERMS 

In this work we are interested in examining the case where a scalar field (j) has a background value (0) — ip{t). 
This field could be, for example, the order parameter in some phase transition or the inflaton field during cosmic 
inflation. For the purposes of these calculations, the scalar field is treated generically. This field is coupled to other 
boson X and fermion ip-)^ fields, which will be called the "catalyst fields" . These fields are in turn coupled to bosonic 
<j and fermionic ijja fields, which will be called the "radiation fields". We will compute the dissipation coefficients 
for the interaction structures listed below. The combination of all the interaction types to be considered encompass 
most of the couplings between scalars and fermions that emerge in particle physics models. In particular, they can 
be combined with appropriately chosen coupling constants to be applicable in Supersymmetric (SUSY) models, such 
as often used in cosmological warm inflation model building [l8| . Alternatively they can be used independently. In 
Appendix VK\ we will give some examples. 

Next we specify the different types of interactions we will consider. The interactions between catalyst and radiation 
fields are: 

1. Interaction between the scalars x and cr: 

£z[x,X^'T,at] = -ftiM[xV2 +x(at)2] -53VxaV ^ (3^^) 

where hi and 53 are coupling constants and M is a mass term (in SUSY models this often is M = g\ipl\f2, 
where ip is the background inflaton field). 

2. Yukawa interaction between the scalar x and fermions t/'o-, V'tr^ 

'C/[X,X^'0<T,'0^T] = -hl\x^i)aPRi'o + Xi^aPLlpa] , (3.2) 

where Pr(l) = (1 i 75)/2 are the chiral projection operators. 

3. Yukawa interaction between the scalar a and fermions ipcr,ipin'4'xT'^x' 

£/[ct, o-''', Vcr, V'o-, V'x: V'x] = -h-i[<y'^i'xPB.i^<y + cripaPLipxl ■ (3-3) 
The interactions between the scalar (f) field and the other fields are: 
1. Interaction between the scalars i/i and x- 

Li[c^,^\x,X^^-9l<l>Ux^X, (3.4) 
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where, without loss of generahty, the background scalar field is associated with the zero mode of the real 
component^ of ^: (f> = {(pi + (p + i(/)2)/-\/2- From Eq. p.4p . the relevant term that will contribute to dissipation 
in the scalar fields effective equation of motion, at one-loop order, is then —g1ip'^\x\'^/2. 

2. Interaction between the scalar (j) and the fermions ip-^, ip^, 

Li[<j>,(j)\ip^,'ip^] = -52[0''''0x-PflV'x + '^V'x-Pi'V'x] , (3.5) 

where, in this case, the relevant term that will contribute to the scalar field's effective equation of motion, is 

-52(¥'/V2)V'xV'x- 

3. Interaction among the three scalars (j), x ^^^d a: 

£/[0,0^x,x^'T,at] = -gi/ii[0t^t^2+^^(^t)2] _ (3 6) 

Note that the interaction (|3.II) is the term dependent on the scalar background field, derived from p.6p once (p 
is written in terms of if and AI is then associated with gxLpl\j2. The interaction term p.6p only contributes to 
the dissipation in the scalar field's effective equation of motion at two-loop order. 

Finally there can be self-interactions for the scalars: 

£7[0,</)t] + Li\x.x^\ + £7[a,at] = A4#t)2 + A^Cxx^)^ + \.{oo'^f . (3.7) 

In the interaction forms considered above, at tree-level </> is only coupled to either the scalar boson x or fermion 
■0^. In the presence of a vacuum expectation value ip^ these can become heavier than the scalar cr and fermion -00-, 
which in turn are only coupled at tree-level with the intermediate fields x. ^-nd V'x ■ -^f the level of model interactions 
we are considering, there are radiative (loop generated) terms that can couple to cr and ^/lo-, which can make them 
also heavy fields. These radiative generated interactions can however be kept small with either suitable choice of 
coupling constants or by imposing a symmetry for the model interactions (like SUSY). These possibilities, together 
with physical examples, are discussed in Appendix [X] In the following sections we will always assume that the 
radiation scalar cr and fermion ip^j fields are lighter than the intermediate scalar boson x or fermion ip^ fields; it is left 
as a model building question what conditions are required to achieve these conditions, but aside from some comments 
in Appendix El those issues are beyond the subject of this paper. 

A. Model Parameters 

The calculation will be done for different thermal regimes depending on the relation between temperature T and 
the field masses m^,m^^,m^,m^^, for the x, V'x' '-'^ ^^id ipa fields respectively: 

• Very low-temperature regime: m^,m^^,ma,m^^ ^ T. 

• Low-temperature regime: m^, m^^ 3> T 3> m^^ ^ip^- 

• High-temperature regime: m^,m^_^,m(j,m^^ <^ T. 

The very low-temperature regime has all field excitations suppressed by Boltzmann exponential factors, and thus so 
are the dissipation coefficients. In the low-temperature regime the interactions produce the two-stage catalyst induced 
dissipation of j7i]. This mechanism leads to dissipation that are not damped by Boltzmann exponential factors, despite 
the fact the x and V'x fields are heavy with respect to T. For the dissipative coefficients, the interactions that give the 
leading T dependent terms have been calculated in [lo| , but here we also present interactions that give less dominant 
terms with respect to T. 



^ Note that here we define the complex scalar fields in terms of their real and imaginary components like for example: <f> = {(j>i + i(f>2)/V2 
and similarly for x and cr. 
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Dissipation coefficients have already been computed for the self-interacting scalar field, in particular in the high- 
temperature regime [Tgl . [20| . In this case, for example in the computation of dissipation coefficients in the local 
approximation, it demands much more stringent restrictions on the model parameters 0, H than dissipation 
coming from (f) interacting with other fields. We will not be interested in the dissipation effects from (f) self-interaction. 
As such we will assume the cf) self-coupling = 0, or we can equivalently assume that is only a classical background 
field which is coupled to the remaining fields (considered as the bath or environment which is coupled to the system 
background field). 

Concerning the magnitude of the many interaction terms considered above, we will restrict to values such that the 
leading contributions to the vacuum polarization terms (from which the decay terms are extracted and, as we will 
see below, enter in the calculations of the Green's functions) will come dominantly from one-loop processes, when the 
decaying field is a heavy field (e.g. either x or ipx)- For instance, we will keep the leading one-loop order vacuum 
polarization terms for the catalyst (intermediate) fields x and coming from the radiation fields a and f/'a- , and that 
leads to decay processes like x — J' f -I- ct, x ~^ "^o- + V'o-j and V'x ~^ + V'o- j with complex conjugate analogous process. 
Decay processes for the light boson and fermion fields can only happen through Landau damping and, therefore 
vanish at zero temperature when on-shell. Decays for the heavy fields involving two-loop terms, like Landau damping 
coming from the fields self-interactions or the bi-quadratic interaction will be considered subdominant with respect 
the one- loop order processes. For this to be valid, the coupling constants in the interactions (13. ip . (|3.2p . p.Sp and 
p.7p should satisfy 



gj < hi,h2 

Ax</ii,/i2- (3.8) 

On the other hand, decay processes involving one light and one heavy field are suppressed compared to the two-loop 
process involving the a self-interaction term with coupling Ao- . The only restriction we impose on the cr self- interaction 
magnitude is it to be small compared with the couplings, which then implies in the constraint 



X,^glhl,hl, (3.9) 

so that temperature corrections to the heavy fields tree-level masses are larger than the temperature correction for 
the cr. This way we can also keep the effective (temperature dependent) masses for x and t/'^ still larger than that for 
the a and ip^ in the high-temperature regime. 

In the low-temperature regime, the above generic perturbative constraints are sufficient for consistent calculations 
at one-loop order with respect to the intermediate fields. In the high-temperature regime, however, computation of the 
dissipation coefficients are known to have issues regarding resummation of higher loop terms [2l| - |23j . We will return 
later on to the issue of resummation when computing the dissipation coefficients in the high-temperature regime. 
There we will present some additional constraints on the model parameters, in addition to (j3.8p . so as to make the 
calculation consistent in that regime. 



IV. DISSIPATION COEFFICIENTS FROM SCALAR AND FERMION INTERACTIONS 

This section presents the specific dissipation coefficients emerging from the interaction terms given in Sect. Imland 
that appear in the (/)-field background effective equation of motion. 

As discussed in Sec. |lTl the effective dynamics for the background scalar field when expressed in a local form is 
described by a dissipative equation of motion of the form Eq. (j2.8p , with dissipation coefficient T expressed in general 
by Eqs. ([21]) and (|2T0)) . 

As mentioned before, we will assume throughout that (f> contributes only through its classical background field 
value. Moreover, these calculations will also be done in the adiabatic approximation Eq. (|2.6p as well as 

^«^"'- (4.1) 

Thus, regarding the interaction terms of the intermediate (catalyst) fields x and V'x presented in Sect. IIIIl we find 
that Y^n, the retarded correlation entering in the definition of the dissipation coefficient T, Eq. (12.91) . becomes 
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where the first term in the above expression conies from the interaction (|3.4p and the second term from (|3.5p . The 
dissipation coefficient T as defined here is derived for a system near thermal equilibrium. As such, the averages in Eq. 
()4.2|) are over a thermal equilibrium distribution. Equation (|4.2|) is represented diagrammatically by Fig. [1] where 
only the relevant contributions to the vacuum polarization (as discussed in the previous section), which enter in the 
computation of the field propagators, are shown. 



Xj A 




FIG. 1: The contribution to the dissipation coefficient at leading order. The thick propagator lines are dressed scalar x (dashed) 
or fermionic i/'x (solid) propagators (shown here including the relevant contributions to vacuum polarization at one-loop order). 

In terms of the interactions given by Eqs. (|3.4p . p.Sp and p.6p . the dissipation coefficient at leading order is [1], 

2 /(72\2 „ ^4 

T = yi^yj j ^^i[Pxi{Po.V? + Px2{Po,V?]nB{po)[l + nB{po)\ 

+ ^ / (^tr[p^^(po,p)']nf(po)[l-nF(po)] , (4.3) 

where and are the spectral functions for the intermediate fields x and defined in Sect. IVII bv Eqs. (|6.ip 
and (|6.3p . and ub and np are the Bose-Einstein (713(0;) — (e*^" — 1)^^) and the Fermi-Dirac (nir(w) = [e^^ + 1)^^) 
distributions, respectively, with j3 = 1/T. 

V. PARTICLE PRODUCTION INTERPRETATION 

The physics involving the appearance of dissipation terms in the background effective equation of motion can 
be given an interpretation in terms of radiation (or particle) production. Such a physical interpretation for the 
dissipation term (12. 9p has been given before 0, [H, [H [13 . We here briefly review this interpretation in terms of 
the field interactions we are considering. Here dissipation can be seen as a result of energy being transferred from 
the (j) field to the radiation bath fields a and ^o- through the excitation of an intermediate field, which in this case 
are either the scalars x or fermions V'x- That this is indeed associated to particle production (ct and/or ■00-) has 
been demonstrated explicitly in In this case, the particle production rate for the radiation fields are only due 
to the interactions of x^'^x with a,ij}a, since the latter do not couple directly to the background field, except by the 
interaction term Eq. (j3.6p . which is however higher order (two- loop), in comparison to the other interactions. The 
particle production rate of the radiation bath particles due to interactions is [2J| 
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fi — Im 



-iu>(p)(t-t') 

2 / dt' '' 

2tj(p) 



-S2i(p,i,t') 



(5.1) 



where w(p) is the particle dispersion relation for the scalars a or fermions ipa and S21 = "^0,21 + 5]^^, 21 is the 
sum of the a and 'i/'o- self-energies (in the Schwinger-Keldshy formalism representation) . n — ha + h^^ is the total 
particle production rate of the radiation bath particles ct, V'cr,^cr from the decay processes x ^ 2fT, x ~^ V'o- + V'o- 
and tp^ ^ a + tp^. The self-energies So-,2i and S^^,2i are obtained (or vice- versa) from the diagrams contributing to 
dissipation shown in Fig. [TJ by cutting one of the internal loop (joining the external) lines for a (in the case of the a 
self-energy) or (for the ipa self-energy). The diagrams representing these self-energies terms are shown in Fig. [21 
An explicit expression for the self-energy is given for example in [23 ]. 



g,9 
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FIG. 2: Diagrams contributing to the self-energy of the radiation bath fields. 
Using the expression of the total radiation energy, 



the dissipation coefficient Eq. (j2.9|) can be expressed as 



T = ^ = — / ——3 [a;<,(p)7i<, -|-a;^„(p)n^J . (5.3) 

This expression directly links the particle production of the bath radiation fields with how the background field 
dissipates its energy, represented by the dissipation coefficient T. 

VI. SPECTRAL FUNCTIONS 

Dissipation coefficients require the spectral functions for the scalar and the fermion fields. For the scalars (where 
j — Xij X2j cTij "'2 denote the real components of the complex scalar fields) we have: 

I X ^ i i 

PAPo, Pj _^2 p2 ^. ^2 ^ iImS,(p) -pI + p2 + m|j . - ilmSj(p) 



4u}j{p)T^{po,p) 



[-pl+ujj{prf + [2c^,(p)r,(po,p)]' 



(6.1) 



where ^^(p)^ = P^ + TTT-j^ j, with j the effective, renormalized mass for the scalar fields, rnj^j = rrij + ReSj. Also 
in Eq. (j6.1|) the decay width, Tj(po, p), is defined in terms of the imaginary part of the self-energy as 
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Note in the adiabatic approximation Eqs. (12.61) and (|4.ip . which under Ues aU calculations here for dissipation coeffi- 
cients, the response time scale is typically associated with the decay widths, r — l/Fj. 
For the fcrmions (where now j = ip-^,ip^) we have that 



^•'^^^ ^ + TTifij + ilmT,j{p) ^ + mflj- — ilmEj (p) ' ^ 

where m^j- is the effective, renormalized mass for the fermion fields, mj^j — rrij + ReSj. The imaginary part of 
the fermion self-energy appearing in Eq. (j6.3p can be written, using rotational invariance, as (1 is the 4x4 identity 
matrix) 



ImSj(p) = 7 po&Q + j.pap + mR.jam'^ , (6.4) 

where the expressions for olq , dtp and are explicitly given in Appendix [Bj 

The imaginary part for the self-energies, appearing in Eqs. (16. ip and (|6.3p are derived in Appendix |B] for complete- 
ness. The results for the decay widths are also given in Appendix |B] For example, for the bosonic interaction Eq. 
(13. ip . it is given by Eq. (IB9p with g\ = 4(/iiM/-\/2)^, while for the Yukawa interaction Eq. p.2p . it is given by Eq. 
(|B12p with coupling factors now defined by g\ = (/i2/-\/2)^, and for the mixed boson-fermion interaction Eq. p.3p . 
it is obtained, likewise, from Eqs. (|B17p - (IB2C)p . with g\ = hi. 



VII. DISSIPATION COEFFICIENTS IN THE LOW AND HIGH TEMPERATURE REGIMES 



An analysis of Eq. (|4.3p indicates that the behavior of the dissipation coefficients in the different temperature and 
interaction regimes is determined by an interplay between the spectral function and the thermal occupation numbers. 
There are two regimes where well defined approximations can be made. One is the regime of fixed T and where the 
couplings of the catalyst fields with the light fields, hi, are going to zero. In this regime, the poles of the spectral 
functions will dominate the integral and this will be referred to as the pole approximation. The other regime is the 
integration region around low energy po and low three-momentum p, and this will be referred to as the low-momentum 
approximation. In the pole approximation, a well defined peak occurs in the spectral function, at an energy around 
the heavy particle dispersion relation. When the temperature is large, the occupation numbers are also large, and 
this approximation is valid. As such, the pole approximation works well in the high-T region. In the low-momentum 
approximation, and at low temperatures, T < m^^m^ , the spectral function still has a well defined peak at around 
the heavy particle dispersion relation, but now this contribution is exponentially suppressed due to the occupation 
numbers in Eq. (14. 3p . However even at low T the occupation numbers in the low energy and momentum region will 
not be exponentially suppressed, and so this approximation dominates over the pole contribution, which is hugely 
suppressed from the exponentially damped occupation numbers. Thus at low T the low-momentum approximation 
dominates. Finally, at very low T when also the decay widths become exponentially damped, both the pole and 
low-momentum approximation can be important, depending on the size of the couplings of the catalyst fields with 
the light fields, hi, and the masses of all the particles involved. 

Both approximations are now examined in detail. In the pole approximation, valid for fixed T and ft. — > 0, the 
general expression for the dissipation coefficient Eq. (14. 3p take on the simplified form. 



2g| f d^p m\^^ 

+ -] (^ r,j.,,,p).^/ "^"^-) - (7.1) 

where Fj is the total decay width, e.g. for xi.2 it is given by the sum of Eqs. (jB9p and (jB12p . while for fermions, e.g. 
F^^ , the fermion ■0^ decay width is defined by the imaginary part of the pole of the fermion spectral function Eq. 
(|6.3p and is given for example by Eq. (IB17p . 
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For the low-momentum approximation, valid for fixed couplings hi and T 
boson fields becomes 



0, the spectral functions for the scalar 



while the spectral function for the fermion fields, Eq. 
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becomes 
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(7.3) 



and then 



tr(p?) ^ 16 



(7.4) 



In this approximation, the dissipation coefficient, Eq. ()4.3p . using Eqs. ()7.2p and ()7.4p . becomes 
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(7.5) 



Note from Eqs. (|7.ip and (|7.5p that the dependence of the dissipation coefficient on the heavy fields catalyst decay 
widths is very different in the two regimes where the pole and low temperature approximations are applicable. While 
in the pole approximation the dissipation coefficient is inversely proportional to Fi, in the low temperature regime 
it is proportional to F|. These different dependencies will be of fundamental relevance when discussing higher order 
contributions to the dissipation as we will analyze later on. Next, we analyze the dissipation coefficient coming from 
the scalar and fermion interactions in the many different energy regimes with respect to the field masses and the 
temperature. 



A. Very Low Temperature Regime 



In the very low-temperature regime, T <C m^, , m^. , m^^, the decay widths for the heavy fields into light 
fields, which can be read from Eqs. (IB9p . (IB12P and (jB17p . with X the heavy field and the decay products Xi,X2 
corresponding to the light fields, are dominated by the zero temperature contributions, 
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(7.6) 



with temperature corrections that are Boltzmann suppressed by the heavy catalyst field masses. 

Both the pole and low-momentum approximation have regions of validity and in both cases, as we will see below, 
the dissipation coefficient associated with a heavy field decay get exponentially suppressed. In general in this very low- 
temperature regime, for small couplings hi, the pole approximation, Eq. ()7.ip holds. The integrals in this expression 
are dominated by low three-momentum with respect to the field masses, to give 



T<m. low h 

Two 



.(1) 







(7.7) 



11 



For large values of couplings hi in this very low-temperature regime, the dominant contributions to the dissipation 
coefficient, Eq. (j4.3p . come from the low-momentum approximation to the spectral functions, Eqs. (I7.2p and (|7.3p . 
with po, IpI ^ "^xi ' "^'Ax ' ™<^i ' ^'^'/'.T ■ Bose-Einstein and Fermi-Dirac distributions in Eq. (|4.3p . the larger 

contribution to the energy integrals now comes from T ~ po- The dissipation coefficient can then be approximated as 



91^ 



■ r2r(2)/r 



(7.8) 



The behavior seen above is directly proportional to the square of the decay rates. If po ^ p, the energy conditions in 
Eqs. (jB9p . (jB12|) . (jB18|) . (jB20p and (jB19|) cannot be satisfied and the the dissipation coefficient vanishes identically, 
T = 0. For pq < p, only the Landau damping terms in the decay rates contribute. But these are exponentially 
(Boltzmann) suppressed, thus the dissipation coefficients in this case are negligibly small. As T — >• 0, T — )• 0, as it is 
expected [l7| . 



B. Low Temperature Regime 



The low-temperature regime is defined as Ti^^ij, m^^^^ <^ T ^ i^R,xii'^R,i! ■ other words, this is the regime 
of low temperature with respect to the heavy catalyst fields but high temperature with respect to the radiation 
fields. In general the masses of the radiation bath fields mcr,m^^ can be neglected. The temperature corrections 
to the effective masses of the heavy fields, due to the light field self-energies, can also be neglected, e.g., ^ — 
+ (4A^ + 93+ h^T"^ /12 « and ~ m^^ -I- /i§T^/12 « "i^^, likewise for the light fermion field ipcn whose 

loop radiative corrections, including the thermal ones, are suppressed by the heavy field masses, even though in this 
case the light fields are in the high-temperature regime. However the light scalar a field effective mass will get a 
temperature correction from self-interaction, m|j ^ ~ -I- XJT^ j'i. 

In this regime one can verify that the dominant contributions to the dissipation, Eq. (j4.3p . come from poj |p| ^ 
mB.,xi ' '^R,il) J which lead to the low-momentum approximation. There is also a regime for fixed T and /i — >■ in 
which the pole approximation is valid. Since the temperature is small with respect the heavy field masses, the three 
momentum integration in Eq. (j7.ip is dominated by |p| ^ mfl^-^. ,mfl_^^ . In the pole approximation the decay 
widths are just given by the zero temperature results in Eqs. (|B9p . (jB12p and (|B17p . The behavior of the dissipation 
coefficient is, however, exponentially suppressed and of the form of Eq. (17.71) . 

For the case of fixed h with T — > 0, the low- momentum approximation is valid, so the spectral functions can be 
approximated as Eqs. (|7.2p and (|7.3p . The dissipation coefficient is now given by Eq. (|7.5p . which then leads in this 
low-temperature regime to the result: 



gWT^^ 1 



where the constants in the above equation arise from evaluation of the momentum and energy integrals in Eq. (|7.5p . 
They have been evaluated numerically to give, Ai ~ 124.9, Bi ~ 129.8, Ci ~ 18.2, Di ~ 4362.6. 

Note from Eq. (|7.9p that at leading order the dissipation coefficient goes as 0(T'^/m^), agreeing with |10jj. The 
dissipation coefficient obtained in this low temperature and low momentum regime is much less suppressed than what 
would be expected from the pole approximation case. In particular, we obtain no exponentially suppressed results as 
in the very low temperature regime. 

The heavy fermion loop contribution is worth noting in Eq. (|7.9p . It is given by the last term on the RHS of that 
expression. The fermion contribution is sub-leading in the temperature in the low-T regime, as compared with the 
term coming from the decay into only bosons. 

This result can be explicitly verified by numerically solving the full expression, Eq. (j4.3p . and comparing with the 
respective approximations in each temperature regime case. This is done in Figs. [3] and SI for choice of coupling 
h = 0.5 (where we take for simplicity hi = h2 = = h in the expressions for the decay widths) as function of 
temperature T, at both low and high temperature respectively. It is clear from Fig. [3] that as T <C m^, the low- 
momentum approximation agrees well with the full expression Eq. (j4.3p . whereas the pole expression agrees poorly 
as expected. The reliability of the results at high temperature, where the pole approximation is applicable, will be 
analyzed in more details in the next subsection. 

From these comparisons, we see in all cases, the low-momentum approximation works well at very low temperatures 
T/to^ ^ 0.04, whereas at higher temperatures T/m^ > 1.0, the pole approximation works well. The breakdown of the 



12 





dissipation T 




0.026gthiip'^T^M'^/m%^^ 


-gl<t''''l>x''x - h2[xH<TPR-4'a + xi^^PL^T] 




^gUHx^X-h^Mix^a^+xi'y^?] 


0.015gf hih^^V^T" My m%^^ 






--^g^Vi'xi'x - hsla^'ip^PRipa + atpxPLlpa] 


0.22gihtTym%,^^ 



TABLE I: A summary of all the dissipation coefficients in the low-temperature regime coming from each of the interaction 
cases involving a heavy intermediate field (all expressions evaluated with m^i = 1^x2) ■ 

low-momentum approximation as T increases arises because the occupation numbers in Eq. (j4.3p no longer become 
exponentially suppressed, and so the contribution from the poles start to dominate. Results for the dissipation 
coefficient, coming from the individual interaction types, in the low-temperature regime, are summarized in Table ID 




0.07 



T/m 

X 

FIG. 3: The behavior of the dissipation coefficient as a function of the temperature (in the low-temperature regime), for gi = (72, 
= /ii = = /i3 = 0.5, rria = O.OOlmx and m^^ = m^, m^^ = nia- 



C. High Temperature Regime 

In this regime, we still have ma,m^^ <C T. As for the heavy catalyst field contributions, their effective thermal 
masses are TOp^ ^. ~ m^. + (4A^ + g| + /i2)r^/12 and TTip^ ^ — rri^ +hfT^/12. Two temperature regions can be useful 

to define, m^.,m^^ <^ T < ^/Vlm^J ^j4X^+g^+l^^,y/\2m^^/hz^ and ^/Vlmj^J g'^ -\- h^, ^/V2m^_^/h^ < T. In 
both cases Eq. (|4.3I) lead to a high-temperature expression for dissipative coefficient, which is dominated by the pole 
approximation and lead to the approximate expression Eq. (j7.1l) . In this temperature regime, leading order one- loop 
expressions for these coefficients can become problematic as we will see below and are only valid in a restricted region 
of parameter space. The reason for this can be seen from how the dissipative coefficient in Eq. (j7.ip depends on 
the field decay widths F^. Since it is inversely proportional to F; and since the F^ are proportional to the square 
of the coupling constants (see the results in Appendix |B]) , loop corrections come with a ratio of the vertex coupling 
squared to the decay width coupling squared. Depending of the magnitude of these two couplings, in some parameter 
regimes higher loop correction terms to the dissipation will be the same order as the one-loop term, thus requiring 
resummation, and in other parameter regimes higher loop contributions will be suppressed. 
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1. Resummation issues 



It is worth recalling at this point the origin of this dependence on the couplings. At high temperature, diagrams 
contributing to the dissipation would have near on-shell singularities if only products of bare propagators were used. 
These singularities are softened by introducing explicit lifetimes for excitations, thus the decay widths F;, through 
dressed propagators. This procedure leads to the result shown in Eq. ()7.1|) with a similar procedure at higher loops. 
The presence of these regulating thermal widths in the denominator then imply entire classes of diagrams, called 
ladder diagrams, which are diagrams with insertions of loops between the external propagators and in the definition 
of the dissipation they can contribute at the same order. This is much like what happens in the definition of viscosity 
coefficients from Kubo formulas as seen for example in Refs. [Ill, [l^, in the case of a single scalar field, and full 
resummation schemes to account for all the corrections occurring at the same order were used. Here we apply the 
outcome of this analysis for our case of field interactions involving light radiation and heavier catalyst fields. 

Two possible examples of topologies of diagrams that can contribute beyond one-loop order are shown in Fig. [5] 
Higher order diagrams can be made with different variations of these topologies and similar ones. 

Using the simple cutting rules explained e.g. in Refs. [HI, H^l; l^t us first consider the contributions in the first 
topology shown by the top diagram in Fig. [5j which are the types of topologies contributing to dissipation due to the 
heavy catalyst scalar field x- Higher order contributions involve adding rungs to the simple leading order contribution, 
which are of the form of the rungs (a)-(e). The analysis from these contributions to dissipation shows propagators 
which can share pinching pole singularities. These are propagators for the fields running on the external part of the 
diagram that can share the same momentum when going on-shell, thus displaying near on-shell pinching poles. As 
mentioned above, these singularities are regulated by the field's decay width, giving as a leading order contribution 
to the dissipation for instance a contribution of order (gf/Fi)", depending the type and number of rungs added to 
the diagram. For example, take for instance the case of a rung of type (a) shown in the top topology in Fig. [S] 
In this case the contribution of the rung (a) produces a multiplicative contribution to the leading order one-loop 
dissipation term that is 0(A^/r^). F^ receives contributions coming from the x self-coupling vertex, ^ O(A^), from 
the coupling with the light a bosons, ~ 0(/if), and from light ^a- fermions, ~ 0(/i2)- this case, this contribution 
can be rendered small compared to the one- loop result for the heavy field provided that /ii, /12 A^. Another possible 
contribution comes from rungs of type (b). Its contribution is 0{g^/T^), which can be rendered small by the same 
considerations taken for type (a). The next case are contributions involving rungs like type (c). Its contribution is 
now of 0{g^/Ta-) The thermal width for a comes from the imaginary part of the self-energy terms for a built with 
the interactions (|3.ip . p.3p . p.6p and its self- interaction term in (|3.7p . The leading order contributions from these 
interactions give F^ ~ 0{Xl) + 0{hl) + 0{hl) + 0{gl) ~ 0{hl,hl). These contributions can also be rendered small, 
compared with the one-loop contribution with external propagators given by the heavy field x, with a suitable choice 
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(a) (b) (c) (d) (e) 




(f) (g) 




FIG. 5: Examples of topologies of diagrams that can contribute to the dissipation beyond the one-loop order. 

of couplings, e.g. 53 ^ /ii,/i2,^3- The next two cases of contributions involves rungs of the form (d) and (e). The 
contributions of the form (d) gives 0{h2/(rx^->P<T))- Since T^^ ~ 0(^1 + ^3) ^^'^ using e.g. the typical parameters 
found in SUSY models displaying the types of interactions we have, the couplings will satisfy Eq. (jA6[) . and we find 
that these contributions can potentially be of 0(1). Likewise, the contribution of type (e) for the couplings shown 
in Eq. (|A6p gives a parametric result of 0{g'^/h'^), which is higher order. The last topologies to be analyzed are the 
ones contributing to the dissipation coming from heavy catalyst fermion field and given by diagrams of the form of 
the bottom diagram shown in Fig. [S] The relevant additional rungs that appear are of the form of the terms (f ) and 
(g) shown in that diagram. The contribution from rungs of the form (f) are parametrically 0(ft.|/(r^^r^_^)). Using 
that ~ 0{h^), and parameters shown in Eq. (|A6p . we find that this contribution is potentially of 0(1). The 
other contribution comes from rungs of the form (g), which gives 0{h^/ (r^^Ta-)) ^ 0(1) using again Eq. (|A6|) . This 
contribution is also potentially of same order as the one-loop result for the dissipation due to heavy fermions in the 
higher temperature regime. The fermion dissipation, however, tend to be numerically smaller than the scalar field 
contributions to dissipation at high temperatures, as we will see from the results to be shown below and from those 
in Fig. ini 

2. One-loop results 

Though in principle all higher loop contributions to the dissipation can be re-summed using a variety of techniques, 
we will not attempt to do this here. The interactions considered in this paper are much more complex than the simple 
case of only a single scalar field, which was considered in previous works performing these resummations, and it would 
require a separate and thorough analysis to examine the necessary resummation. In any event, it is important to 
recognize that in contrast to the original analysis made e.g. in Refs. [lll,[12| that had only a single coupling parameter, 
in our case where we have multiple coupling parameters and, as discussed above, we now have the added freedom to 
choose certain parameter regimes where some of higher loop contributions can be suppressed and no resummation 
for them are in principle necessary. This is true for example when computing particular contributions of fields, where 
for example, for those classes of diagrams that involve intermediate heavy field propagators and exemplified in the 
previous subsection. The cases that need special attention with respect to resummation are those where light field 
components run in the external part of the diagrams contributing to dissipation. This is e.g. the cases involving the 
inclusion of rungs of type (c) and (e) in the top topology shown in Fig. [SJ or the rungs (f) and (g) in the second 
topology of the same figure, and seen in the parametric analysis performed previously. 

Below we show some results for the dissipation at one-loop order at high temperatures, but these results are only 
meant to be phenomenological ones, providing an estimate for the dissipation, since they do not account for all the 
possible higher order loop contributions that they may receive. 

In Fig. Ulare shown the result for the total dissipation coefficient, for the case of summing all the contributions, 
for temperatures up to two hundred times the intermediate field mass. Note that the low-momentum approximation 
quickly overestimates the true result as the temperature increases. The pole approximation becomes very good already 
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at T ~ my^,m^^, but tends to slightly lose accuracy at higher temperatures, where the simple pole approximation, 
given in terms of the particles dispersion relation, becomes less accurate. This is more critical in the case when 
only decay into bosons are considered (e.g. with decay width given by Eq. (jB9| ). and the temperature is such that 

T^^^ /m^ > 1 (note that this cannot happen in the case of decay into fermions, where for large temperatures the decay 
width tends to zero (on-shell) due to Pauh blocking [25|). 

It is also worth showing the individual contributions to the dissipation coming from the different intermediate 
catalyst fields and the different decay channels available to them. This is shown in Fig. [51 




Note from Fig. [5] that for the cases of interaction of a scalar intermediate field with light scalars and interaction of 
an intermediate fermion with light scalar and fermion, the high temperature behavior is linear in temperature. Also, 
the dissipation from the interaction of an intermediate scalar field with light fermions dominates in the regime of 
temperature where T < m~^^/h2, but for T ^ m^/ft.2i the interaction of a intermediate scalar with light scalars 

dominates. For T ^ m-^jh^ the contribution to dissipation coming from the interaction of an intermediate scalar 
field with light fermions tend quickly to become subdominant with respect to both the interaction of the intermediate 
heavy scalar x with the light scalars a and the intermediate heavy fermion with the light scalars and fermions. 

Approximate semi-analytical fitting formulas for the dissipation, valid in both high-temperature regimes, with 
respect to the masses of the intermediate fields, m <^ T < m/h and T ^ m/h, and expressed in terms of the 
thermal masses for the intermediate scalar and fermion fields, can be obtained by considering separately the different 
interaction forms. Again, we cautioned that these results at high temperature are only meant to be interpreted 
phenomenologically, sincethey are only derived at leading one- loop order. The different cases are listed below (for 
simplicity all expressions below were evaluated with m^-^^ — m-^^): 

a) Interaction of a heavy scalar into light scalar bosons, Eq. (13. ip : 
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where the thermal mass (T) is here defined as (T) ci 
b) Interaction of a heavy scalar into light fermions, Eq. 
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g^T / 12, due to the interactions shown in Eq. (13.1 
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where the thermal mass 'm^{T) is here defined as m'^{T) ~ + h^T'^ /12^ due to the interactions shown in Eq. (13.21) . 
c) Interaction of a heavy fermion into a hght scalar boson and light fermion, Eq. p.3p : 
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(7.12) 



XT) 



due to the interactions shown in 



where the thermal mass m^p^ (T) is here defined as rn^ 
Eq. (ig. 

As noted above from the results of Fig. [51 the case of interaction of the intermediate scalar field x with boson 
radiation fields a leads to growth linear in temperature of the dissipation coefficient at high temperature. This differs 
from the case of the scalar x self-interacting, which found 1/T behavior [l3|, [l3, E^- The reason is, for the self- 
interacting case, the decay width at lowest order was at two-loops, whereas for the above case involving a separate 
radiation field, the decay width at lowest order is at one-loop. At one-loop the decay width is proportional to an 
imaginary part of the self-energy that is 0(/i^to^), while at too-loop order, there are contributions [ij,[lll of 0{g^T'^), 
coming from the interaction .93X^0-^ and O(A^r^), coming from the x field self-interaction, X^x^- This difference in the 
decay rates at one-loop order and two-loops, results in the differing temperature behaviors observed for the dissipation. 
This difference in the leading behavior of the decay rates at one-loop and two-loop orders also imply that the one-loop 
contributions are dominant provided the temperature is not too high, that is, T < min{him^/g^,him^/X^), which 
gives an additional constraint for the consistency of the derivations done here based on the contributions of one-loop 
polarizations only. 

A crosscheck that the above estimates are at least parametrically correct even in the absence of a full resummation 
treatment can be made by noting the close relation between the general dissipation coeflacient expression, Eqs. (|2.9p . 
(|2.10p . with those of the viscosity coefficients derived through the Kubo formulas [21I, One observes that aside 
for differences in powers of momentum in the defining integrals, the expressions are similar. This similarity can be 
used to relate the dissipation coefficient with the known results for the viscosities, whose results, found for example 
in Refs. (2ll . [2^ . account for the full resummation of ladder diagrams, just like the ones also contributing to the 
dissipation coefficient at high temperatures. This type of approach was used for instance by the author in Ref. (26j 
to relate the moduli decay rate to the bulk viscosity coefficient. This was possible because of the interaction form 
used, which involved a derivative. This allowed turning the calculation of the dissipation coefficient in that reference 
to the analogous computation of the bulk viscosity through use of the Kubo formula. In our case, following the 
lead of Ref. 26] , the proper analogy would be instead to relate our expressions for the dissipation with that for the 
shear viscosity (and properly accounting for the difference in powers of momentum in the integrals). This analogy 
was in fact raised before in Ref. [l3| . In the simplest case of the interactions of the catalyst scalar field x with only 
the radiation scalar field a, this allow us to obtain an estimate for the shear viscosity as ^shcar ~ /{h\m?,). The 



different temperature dependence with respect to the standard scalar field result obtained in the theory 21 
Cshear ^ / , comcs from the different contribution to the decay width F in each case (as we noted in the previous 
paragraph, in the Xcf)'^ model F comes from the imaginary part of the two- loop self-energy term, which is ^ A^T^, while 
for the scalar catalyst-radiation interaction we consider here, it is dominated by the one-loop self-energy computed 
in the Appendix B). Relating the shear viscosity expression with our expression for the dissipation coefficient in this 
case (and accounting for the differences in factors of coupling constants and powers of temperature), in the high 
temperature regime we obtain T ~ gfT/hi- This is parametrically the result we found in Eq. (|7.10p . Similarly, 
we could also relate the dissipation coefficient obtained for the other interaction forms. But this is still only meant 
to provide a parametric estimate for the dissipation coefficients, since they are not the complete resummation of 
all contributing higher order terms. Thus, we refrain ourselves from following any further this similarity between 
dissipation and viscosity coefficients in our case here, because it would also require a full calculation of the viscosity 
coefficients as well, which is beyond the scope of the present paper. 



3. Contrasting the derivation of the high versus low temperature dissipation coefficients 

At this stage it is useful to contrast the derivation of the dissipation coefficients at high temperatures from this 
subsection with those obtained in the low temperature regime of subsection I VII Bl One recalls from the calculations 
done in the previous subsection IVII B[ that in the low-temperature case, where what we called the low-momentum 
approximation is valid, the dissipation coefficients are proportional to the square of the decay widths. This is evident 
from the behavior of the spectral functions, such as in Eqs. (|7.2p . (|7.3p and from the explicit expression for the 
dissipation coefficient in that regime, Eq. (j7.5p . This very different behavior of the dissipation coefficients in the low- 
temperature regime is why they do not suffer from the problems seen for the high temperature regime. Due to this 
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fact, it also means no resummation is required in the low-temperature case. Figure |3] provided one numerical example 
of parameter and low-temperature regimes where this low-momentum approximation is realized. As no resummation 
is needed, the one-loop expression for the dissipation coefficients in this regime alone provides reliable results. 

It is also instructive and easy to understand the low-temperature dissipation coefficients by looking at the diagrams 
shown in Fig. [5] The dissipation coefficients always carry catalyst propagators in their upmost extreme points. Any 
higher order contribution to the dissipation coefficient would then have at least four external heavy catalyst field 
propagators. In this regime, the catalyst fields are at low-temperature, but the light radiation fields are at high 
temperatures. Thus, the rungs made with the light fields will be the ones that will have the most infrared sensitivity 
and that can produce 0(1) contributions as seen from the analysis made following Fig. [5j But accounting for the at 
least four external heavy catalyst fields in these diagrams, from the highest order contribution from them, we are led 
to contributions of order 0(rj) x 0(1), which are still of subleading order compared to the leading one- loop results 
for the low-temperature regime given in Table U which recall are O(r^). In the other case, two of these propagators 
can be just Feynman propagators, which will lead to an overall suppression factor 0(T'*/m^) at low temperatures. 
Since the dissipation coefficients in the low-temperature regime have direct application to the warm inflation scenario, 
where they were shown to be the most useful coefficients in model building calculations , our results therefore 

have immediate utility. 



VIII. APPLICATIONS 



The dissipation coefficients computed in the previous sections have numerous applications in cosmology and particle 
physics. All these coefficients have been computed with the system near thermal equilibrium and under the adiabatic 
approximation Eqs. (|2.6p and (|4.ip with the microscopic timescale r in these equations identified with the inverse 
decay widths F^. In the adiabatic approximation all microscopic physics, that governs the creation of dissipation, 
operates much faster than any macroscopic changes that the dissipation affects. In particular, the consistency of the 
adiabatic approximation requires that the decay widths that enter the dissipation coefficients must be much bigger 
than the motion of in the background scalar field and changes in temperature [i^ . 

r.»|,f (8.1) 

In the cosmological setting, in addition the adiabatic approximation requires 

r. » H, (8.2) 

where H is the Hubble parameters. When the conditions (|8.ip and (|8.2p apply, we can readily use the results for the 
dissipation we have obtained. Below some systems are discussed, where these conditions apply and dissipation effects 
are important. 



A. Thermalization phase ending reheating 

In the standard inflation picture [l|, during inflation the universe supercools. In order to end inflation and enter 
the subsequent radiation dominated regime, a reheating phase is required in which coupling of the inflaton to other 
fields becomes important and coherent oscillations of the inflaton induce particle production in these coupled flelds 
^. There has been exhaustive study of the phases of reheating. In many models there is a flrst phase, preheating 
27|, where there are large oscillations of the inflaton, which through parametric resonance lead to explosive particle 
production. Once the inflaton oscillations decrease, this leads to the second phase of small coherent oscillation leading 
to a perturbative reheating phase In all reheating processes, there is one flnal phase, where the inflaton oscillation 
motion ends, the radiation that has been produced is thermalized, and then the radiation dominated epoch begins. 
This can be called the thermalization phase of reheating. In a generic reheating situation the oscillations of (j) are 
damped due to the Hubble damping term, thus generally is a decreasing function of time. Thus at some point 
the adiabatic condition Eq. (18.11) is realized for particle type i. And depending on the strength of the coupling, the 
other adiabatic conditions Eq. (j8.2p also are realized. This is the adiabatic regime where the macroscopic motion of 
the inflaton field and expansion of the universe becomes slow relative to the microscopic motion. Once that occurs, 
and provided (f> is coupled to other flelds, its evolution equation has the form 



+ (3i/ + T),^ + — = 0, 



(8.3) 
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where depending on how (f> is coupled to the other fields, T are some combination of the dissipative coefficients 
computed in the earlier sections. Here the dissipative term will lead to particle production, with the radiation energy 
density governed by 

Pr + 'iHpr^T^'^. (8.4) 

Thus, this equation will control the final temperature of the universe after reheating. 

The dissipative coefficient already requires the presence of a thermal bath during reheating. In general this is 
not difficult to achieve, and preheating for example can serve to that end. Being a non-perturbative process, even 
particles heavier than the inflaton may be produced, which then decay into lighter degrees of freedom. After these 
particles thermalize, the initial thermal bath is produced. However preheating does not always occur and even when 
it does, typically it is not enough to convert all the vacuum energy from inflation into radiation. The next phase is 
standard reheating, in which the perturbative decay of the inflaton occurs. This phase is generic and produces an 
initial thermal bath. Eventually as the inflaton motion slows and the thermal bath grows the flnal thermalization 
phase starts, in which particle production occurs through dissipation governed by Eqs. (|8.3|) and (|8.4[) . 

That the adiabatic condition Eq. (18. ip is reached at some point during the oscillations is practically ensured in most 
of the models of inflation. For example, for the chaotic inflation model A(^^, during the oscillations it is found that 
(j) slows as ~ and ~ a~^, where a{t) is the cosmic scale factor |28j . Thus (j)/(t) ~ is decreasing, and 
provided that at least F; sa const, or growing, which is often the case, then eventually the adiabatic condition will 
be satisfied. On the other hand if the energy density of the oscillating inflaton behaves as matter, like in quadratic 
chaotic models or in hybrid inflation, then (p/cj) ^ const.. Nevertheless, generically one has (j}/(f> ~ m^, being the 
inflaton mass and the typical value for the frequency of the oscillations, and the condition to start in the adiabatic 
regime reads F^ > m^. This is not difficult to fulfill for the decay rate of a heavier field than the inflaton, with for 
example F^ ~ /i|mi/(87r) being the standard perturbative decay into massless fermions for a particle type i with mass 

In other inflationary scenarios, the inflaton does not undergo a series of damped oscillations about a potential 
minima, but inflation is followed instead by a kination period [20|, and reheating takes place through gravitational 
particle production [30| . If kination lasts long enough, at some point radiation takes over the kinetic energy density of 

the inflaton and the universe will enter the radiation dominated epoch with a temperature Tjnj ~ IQpNs^'^ GeV, where 
Ns is the number of scalar degrees of freedom gravitationally produced. However, the same gravitational mechanism 
also produces a stochastic background of gravitational waves, and the overproduction of gravity waves is one of the 
potential problems of this kind of scenarios. The Big Bang Nucleosynthesis (BBN) bound on the energy density of 
these gravity waves translates then into a lower bound for Ng > O(IOO). On the other hand, during kination the 
inflaton fleld only evolves logarithmically with the scale factor, and then ij)/(p ^ a~^. That is, the adiabatic regime 
can be quickly achieved. Kination can end before the gravitationally produced radiation takes over, and it will be 
followed instead by radiation production through dissipation. The shorter kination period means less gravity waves, 
that could be kept within BBN bounds without any further constraint on the parameters. 

B. Warm inflation 

The standard picture of inflation is where the universe supercools during the inflation phase, since there is no 
particle production. Such a cold inflation picture assumes the inflaton has negligible coupling to other fields during 
the inflation phase. However in a realistic particle physics model, the inflaton fleld could be coupled to other fields. In 
this case an alternative dynamical realization of inflation is possible, which is warm inflation, where particle production 
occurs concurrent to inflationary expansion 0, [sil, H^l , with particle production occurring due to the interaction of 
the inflaton with other fields. The presence of particles then leads to a sustained heat bath. This heat bath influences 
the dynamics of dissipation as well as acts on the inflaton fluctuations. The statistical state of the heat bath could be 
any nonequilibrium state, and a fleld theoretic dynamical calculation is necessary. However in treatments up to now, 
an assumption of thermalization is used, which is then checked self-consistently in the calculation. In such a case, 
the fluctuations of the inflaton, which are the primordial seeds of density perturbations, are thermal in origin {34^- 
Several models of warm inflation under this picture have been studied 6-9, 15, 16, 18, 33, 34]. 

Since inflationary expansion requires the (j) motion to be slow, the adiabatic approximation Eq. (18.11) can be 
applicable depending on the coupling of (j) to other fields. In such a regime, and also provided the second adiabatic 
condition Eq. (|8.2p is satisfied, the evolution equation for the inflaton is governed by Eq. (|8.3p . which in the slow-roll 
regime becomes 



(3i? + T),^ + — ~ 0, 



(8.5) 
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with particle production governed by Eq. (|8.4p . Under these conditions, along with the General Relativity condition 
p„ > Pr, where is the vacuum energy density, a warm inflationary expansion occurs. Moreover typically m^, 
m^^ 3> -ff, in which case the dissipative coefficients computed in the earlier sections can be applied to cosmology. 
Under these conditions and depending how <j) couples to the other fields in the underlying particle physics model, T 
in the above equation represents some combination of dissipative coefficients computed in the previous sections. 

The most successful models of warm inflation [S\ use the two-stage mechanism [7|, in which the inflaton field 
interacts with heavy intermediate fields, e.g. like through the interaction terms (j3.4p and p.Sp . These intermediate 
heavy fields, which can be either bosons or fermions, in turn interact with light fields, which can also be either bosons 
or fermions, with interaction terms realizing the two-stage mechanism 0] such as Eqs. (13. ip . p.2p or p.Sp . At leading 
one- loop order, these decay rates can be expressed e.g. in terms of Eqs. (|6.2p and (|6.4I) . for the interaction terms 
(IXTD . (IX^ and 

Note also that since in warm infiation the heavy intermediate fields have masses m^,m^^ 3> H and even the light 
fields, cr, will generally have masses larger than H, the effects of the expansion can be neglected in the expressions 
involving the respective field propagators and, in particular, the expressions for the dissipation coefficients can be well 
approximated by their Minkowski expressions. 

Once the distinguishing feature of the dynamics in warm inflation is due to dissipation, knowing precisely the 
dissipation coefficients, their temperature dependence and their regime of validity is fundamental. This was recently 
exemplified in the Refs. [ssl l36j . where it was shown how the temperature dependence of the dissipation coefficient 
can affect the power spectrum of perturbations. The results for the dissipation coefficients and the derivation of 
their explicit temperature dependence as we have obtained here becomes then of fundamental relevance in the precise 
determination of the power spectrum in warm inflation. 

Dissipative effects might also be relevant immediately after warm inflation ends, during the flrst stages of the 
reheating process, when different production mechanisms are at work. For example a network of cosmic strings is 
formed at the end of the so-called D-term SUSY hybrid inflation [s^], in which a U{1) symmetry is spontaneously 
broken. If cosmic strings form after a period of warm inflation, the transition will take place with the fields place in 
a thermal bath, and/or with dissipation still operative. This may affect the rate of formation of the cosmic strings 
network, and their amplitude contribution to the CMB. Observations limit their contribution to be less than a 10 % of 
the total amplitude [38| , which translates as usual into constraints on the model parameters. Typically the constraint 
comes from having an inflation scale that is too large, which translates into too large of an amplitude of the primordial 
spectrum generated by the cosmic strings. The constraint can therefore be relaxed by a change of the potential which 
leads to a decrease of the inflationary scale. For example by considering a SUSY model with a non-minimal Kahler 
potential as in Ref. (39l|. Similarly, given that warm inflation typically requires a lower inflationary scale, it may also 
help in reducing the contribution to the spectrum of the cosmic strings. 

If the symmetry broken during the phase transition is instead a discrete symmetry, this will lead to the formation of 
a network of Domain Walls (DW). When the DW are stable, these pose a severe cosmological problem, as the energy 
of the domain wall network will tend to dominate the total energy density, and they induce large T fluctuations in 
the CMB [Zoj. These can be avoided if at the time of the transition one of the vacua is favored oyer the other, for 
example by having a tilt in the potential, or a biased initial configuration, towards one of the vacua 41] such the DW 
evolve and annihilate after formation. Other solution can be to dilute them by having a extra period of inflation after 
the phase transition; of the order of 0(10) e- folds of rapid expansion would be enough to dilute whatever unwanted 
relics have been produced after a normal period of inflation. This can be achieved for example with thermal inflation 
[4^ , where inflation is driven by the T corrections restoring the symmetry and keeping the field in the false vacuum at 
a constant vacuum energy. Similarly, given the appropriate interactions in the model, once the adiabatic conditions 
are fulfilled for having inflation driven by dissipation, it may be enough to have some e-folds of warm inflation to 
dilute the network of DW. There will be less restrictions on the parameters values given that one does not require a 
full period of inflation and does not demand to comply with the CMB constraint on the primordial spectrum. Notice 
however that while thermal inflation is driven by the T corrections due to direct couplings of the inflaton to light 
flelds, in warm inflation the extra bout of inflation will be driven by dissipative effects in a thermal bath, while the 
fleld is moving towards the true minimum of the potential, but with the thermal corrections to the inflaton potential 
being negligible. 

C. Second order phase transitions 

Several relativistic quantum field theory systems exhibit second order phase transitions. U(l) model of scalar 
field coupled to gauge field is a common example (43j . The electroweak phase transition in certain regimes, when 

iTiHiggs ^ n^gauge has been argucd exhibits a second order transition for many models [44| . The chiral phase transition 
is also believed to be second order in QCD for two flavors of massless quarks [1^. Several other examples of second 
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order phase transitions in quantum field theory can also be found [4^ . 

The study of second order phase transitions often focus on the phenomenological Landau- Ginzburg mean field 
theory such as in [i^] or utilize a Langevin equation as in [i^ . Most of these treatments rely on symmetry arguments 
and general principles to analyze the transition, but fall short of a first principles calculation. For that the full 
treatment of all interactions is required, which would then generate dissipative effects of the scalar order parameter. 
For that, the coefficients computed in this paper are applicable. 

IX. SUMMARY AND CONCLUSIONS 

In this paper dissipative coefficients have been computed in the near thermal equilibrium and adiabatic approxima- 
tion, Eqs. (|8.1|) . for all combinations of interactions between scalar and fermion fields. There are several new results 
here. In particular, several new interactions presented in Sect. IIIII have been considered, which have not previously 
been treated in the literature for these types of calculations. For the dissipation coefficients at low temperature this 
has led to sub- leading terms in temperature in Eq. (j7.9l) . At high temperature, the interactions structures we have 
studied lead to linearly growing dissipation coefficient, which differs from the 1/T behavior found in earlier papers 
[3 [3 [3; where only scalar field self-interaction and interaction with another light boson field were considered. 

At high temperature, as is well known and we have noted, there are resummation issues with the derivation of the 
dissipation coefficients at leading one-loop order. In particular there are higher order loop terms giving contributions 
that are of the same order as the one-loop term. We have discussed that these higher loop terms can in principle be 
treated in the same way as in the derivation of viscosity coefficients, which also show the same need for resummations 
and have been treated in the literature. However, due to the complicated nature of the multifield interactions we have 
used in this paper, this is expected to be a much more difficult task. Even so, we have provided estimates for the 
dissipation coefficients in the high temperature regime, which should be interpreted only phenomenologically, giving 
an order of magnitude estimate for them. 

On the other hand, in the low-temperature regime, where the catalyst fields have masses much bigger than the 
temperature scale while the radiation fields are at high temperatures, we have shown that the dissipation coefficients 
in that regime do not suffer from the same resummation issues and the results we have provided in this regime are 
precise. This is a direct consequence of multi-field type interactions treated here. The freedom in coupling and mass 
parameters avails itself so that in certain parameter regimes, calculations only up to one-loop order of the coefficients 
are adequate, without requiring resummation of higher order terms, as usually required in the evaluation of these 
coefficients in single scalar field self-interaction or gauge field cases. In particular by working with multiple types of 
interactions involving intermediate fields and light fields, two independent coupling parameters occur, thus allowing 
for regimes where higher order contributions are sub-leading. This is a huge advantage in the use of these type of 
interactions as compared for example with computing transport coefficients in the single field case. 

There are several applications of these coefficients, and the interaction structures we have studied in cosmology, as 
explained in the last section, but other applications can also be found in the context of heavy ion collision, relativistic 
phase transition and any other case that necessarily involves dissipative processes related to a background or order 
parameter field used to study the dynamics. Many treatments using these coefficients can be found in the literature 
for study of warm infiation, chiral condensate in heavy ion collision, and study of phase transitions. In Sect. I Villi 
we have identified a few more new applications. In particular, we have observed that the final stages of reheating 
after inflation generically will lead to an adiabatic evolution of the inflaton field. During this period, which we have 
termed the thermalization phase after reheating, the dissipation coefficients computed here and related past papers 
are important, and will assist the inflaton in ending its oscillations, exiting entirely from reheating, and flnally leaving 
the universe in a thermal equilibrium state. In a different direction, for infiation scenarios ending with a kination 
period with reheating through gravitational particle production, we have identified that dissipation effects could be 
important, thus decreasing the length of this period, which can have consequences on gravity wave production. 

We have computed the contributions to dissipation from various types of interaction, involving fields with very 
different mass scales. First, since the sorts of interactions we consider in this paper are common in many models, 
the contributions from them to dissipation would still need to be computed to understand how significant they are, 
and this has been achieved here. Moreover, there are very few exact limits or approximations in which dissipation 
coefficients can be computed. Table U contains results from this new low-momentum approximation, first developed in 
[T^ . As such, the results in Table U might prove useful as checks for numerical algorithms for computing dissipation 
coefficients. 
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Appendix A: Model building 

In this Appendix we give some model examples for the interactions considered in section IIIII We have focused 
on a scenario with several scalar fields and fermions; one of the scalars, 4>, acquires a background value ip{t), and 
couples to a set of boson x and fermion tp.^ fields. The latter acquire masses through their interaction with (j). The 
X and fields in turn couple to another set of boson a and fermion tp^ fields, which are kept as light fields. By 
definition therefore this set of light fields cannot couple to 0, otherwise they will pick up a mass similarly to x and 
^p■^. The coupling between the scalar field with non- vanishing background value and a set of light fermions can be 
easily forbidden by imposing an U(l) global symmetry. On the other hand, the coupling between the (p field and the 
light scalars can be avoided in a supersymmetric model plus a global or discrete symmetry. Supersymmetry protects 
the light sector from getting large mass corrections. 



1. Supersymmetric model 

A supersymmetric model with the pattern of interactions considered in section IIIII is given for example by the 
superpotential: 

W = g^X^ + hXY^ , (Al) 

where $, X, Y are superfields with scalar and fermion components given by (0, ip^,), {Xy V'x) (cr, ipa) respectively. 
The scalar potential is given by: 

V = g^lxf + h'\a\^ + VlxH^P + ^ghRelcP^x^a^] + Ah'\x\'W\' , (A2) 
whereas the Yukawa interactions are given by: 

- £f = 2g(p-4}-^PL'4>x + 9X'4>cI>Pl'iPx + '^hx'ipcPLipa + ha^p^PLipa + h.c. . (A3) 
Comparing with the interactions given in section lllll we have for the couplings: 

51 - 32 = 2.9 , (A4) 
53 = 2hi ^h2= 2/13 = 2h , (A5) 
Ax = g^ Xa^h'^y = , (A6) 

and once (p gets a background value (p/\/2, the mass parameters are given by: 

M = V2gLp , = m^^ = m^^ = V2gLp , (A7) 



with nia = = m^^ = m^^ = 0. We recover the interaction terms given in section IIIII plus one additional 
interaction between the boson one heavy fermion ip-^ and one light fermion ip^. Due to the presence of the massive 
fermion in this term, its contribution to the dissipation coefficients will be subdominant (suppressed by extra powers 
of m^) with respect to that given by the interaction of x with two massless fermions ip^. 

The conditions on the couplings given in Eq. p.8p are fulfilled by choosing h 1/2 and <C h. On the other hand, 
the quartic coupling is the same than = /12 = Therefore, in the high-T regime, the temperature dependent 
corrections to the masses of a and x are comparable. On the other hand, in the low-T regime with m^/T > 0(10), 
keeping the a field light when including thermal corrections due to its self-coupling only requires = m'^/T'^, 
which is fulfilled with /i < 1/2. 

So far we have assumed that SUSY is unbroken except by thermal corrections; i.e, from the Lagrangian we have 
M = — TOj/)^, and mg- = 0. However, thermal corrections will split the boson-fermion masses for the X superfield, 
with 

Ami - l-lx - -l^J = '-^T' , (A8) 



and this is in turn can induce a mass term for the light field a at one loop. 
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Therefore, an additional condition on the couplings has to be imposed for the cr field to remain light, (5m^ <C 
the low-T regime: 

' « 1 . (AlO) 



-.2 



3(4^)2 V 9 

Taking into account that we have already required h? 3> g, this reduces the allowed range for the couplings, but still 
with a mild hierarchy of the order of h?' /g ^ 0(10) all the constraints are fulfilled. On the other hand, SUSY can be 
already broken independently of thermal corrections, and in this case we also need to impose Am^/T^ ^ {Att)^ g^ / . 
In addition, SUSY breaking may induced additional trilinear couplings in the scalar potential such that now M > m^. 

To summarize so far, in a SUSY model where the interactions between heavy and light fermions are given in terms 
of two coupling g and h, the consistency conditions for the expressions computed in the low-T regime given in Table 
Uto be valid are: 

/i<l/2 , — ~O(10). (All) 
9 



2. Non-supersymmetric model 

In a non-supersymmetric model, scalar masses are not protected, and with the couplings and IxPl'^'P 

the potential, a mass term for a may be radiatively generated through 

an effective vertex |crp|(/)p. Thus, although symmetry breaking patterns can be chosen such that a sector of the 
boson fields remains massless (the Goldstone bosons), those will couple to the cj) field. 

Instead, the scenario worked out in this paper will correspond to a model with two coupled scalar fields, (j) a-nd Xj 
and two sets of fermions: ipx couples to (j) a-nd it is massive, whereas ipo- only couples to x- Yukawa interactions are 
therefore: 

- -^.F = 520V'x-PlV'x + h2x4^c,PLil^c, + h.c. , (A12) 

and the scalar potential is given by: 

V = XM' + KWt + .9?lxn0P + + /4lxl' • (A13) 

Taking the mass parameters /i^ < and /i^ > 0, the field <j) acquires a vacuum expectation value ip/\/2, and the 
scalar fields get field dependent masses = iJ,^+ gl^p^ and = fJ,'^ + 3X^(f^. In addition, we impose ^ m^, 
i.e., Xif, <^ gl- Therefore, the field is a light field with a small self-interaction. In the absence of any other light field, 
dissipation is then given by the fermionic contribution. In the low-T regime this is the second line in Table H] and in 
the high-T regime, it is given by Eq. (j7.1ip . 



Appendix B: Decay widths 



Below we give some of the details of the calculation of the relevant decay widths used in the expressions for the 
dissipation coefficients. The particle decay widths can be expressed as usual in terms of the imaginary parts of the 
self-energies diagrams [ll| . The relevant self-energy terms in our calculations come from one- loop processes and these 
can all be analytically derived in detail. There are three main processes we consider: (1) decay of a scalar boson to 
scalar bosons; (2) decay of a scalar boson to fermions and (3) decay of a fermion to a scalar boson and a fermion. 



1. decay of a scalar boson to scalar bosons 



Let us first consider the decay of a scalar boson X to two other scalar bosons, e.g. X ^ Xi +X2, with an interaction 
verte x g iven by gxXXiX2. The imaginary part of the self-energy for the leading order contribution is given by (see 
also [13) 

ImEx(po,p) = (e'^^° - 1) J j^;^nBipo - ko)px^{po - ko,p -k)nBiko)px2{ko,'k-) , (Bl) 
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where ub is the Bose-Einstein distribution function and pXi{kQ,]i.) is the free spectral function for the Xi (i=l,2) 
scalar field, 



(B2) 



with a;Xi(k) = y + • "^^^ Dirac-delta functions in Eq. (IBip give the different implicit absorption and decay 
processes involving the X, Xi and X2 particles, including the different branch cuts. Taking for example po > 
(for pq < there is a change of sign in the expressions) and also assuming > "nix^ -i there are two processes 

contributing: one satisfying pq = (p ~ k) + 0^X2 (k), for Pq — > {mx^ + 171x2)'^, corresponding to decay 
of an off-shell X particle into on-shell Xi and X2 particles, and the other for po + wxaCk) — utxiip — k), for 
— |m^^ — mj^J < Pq — < {rrixi — mx2)'^, corresponding to Landau damping, where an off-shell particle scatters 
with two on-shcll particles from the heat bath. Explicitly, for the decay contribution 



ImE^<=-y(po,p)=5i7r J 



(2^)3 4a;xi(p-k)u;x2(k) 
and for the Landau damping contributions, 



[H-nB(a;xi(p- k)) +nB(wx2(k))] (5(po-wxi (p-k)-a;x2(k)) , 

(B3) 



d^k 



1 



(2^)3 4a.Xi(p-kVx2(k; 



[-715(^X1 (p - k)) -|-ns(wx2(k))]<5(po - wxi(p - k) -^^^^^(k)) . 

(B4) 



The angular integration in Eqs. (|B3p and (IB4p can be done using the Dirac-delta functions. Using the formula. 



S{x - Xi) 



\dL\ 
I dx I 



f{x^) = , 



(B5) 



we obtain for example that 



S{po - uj„(k) ~ uj^{p - k)) = I ^'' ^(cosg - cos 6*0) , 

|p||k| 



(B6) 



where cos 6*0 is the solution of po ~ yP^ + k^ — 2|p||k| cos^o 
momentum integration to the values such that — 1 < cosSq 5: 1, or fc_ < < where. 



., f I , , / m\r — m 

± 4 r V -pI + p2 



9 ? \ 2 

nix-^ — 1^X2 
-pI + p2 



4my 
-pI + p2 



-,1/2^ 



(B7) 



and the integration over the momentum in Eqs. (jB3P and (jB4p is now trivial. Defining u}± by 



(B8) 



and by the definition of the decay width, Eq. (|6.2I) . the (one- loop) result coming from the X ^ Xi + X2 process 
involving scalar bosons is 



^x{po,p) = gx 



1 



327rwx(p) 



— h ^In 



1 — e 



1 - 



P0-" + 



5i 



T 



327rwx(p) |p 



■In 



1 — e T 
1 — 



, , .„ _ . . Oipl - p^ - (mx, + m,x2f) 

9{-pl + p^ + imx,-mx2f) , (B9) 



1 - e- 



1-e 



_ P0 + " + 



where has mass dimension two. 
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2. Decay of a scalar boson to fermions 

The second process considered is the decay of a scalar boson X to two fermions. We follow for simplicity the same 
notation as used in the previous case, but now Xi and X2 as a fermion-antifermion pair, which has a Yukawa coupling 
gx with the scalar boson X . The imaginary part of the one-loop X field self-energy is now given by 



^^^nF(po - fco)ni=^(fco)tr [pxAva - fco, P - k)^^^ (fco, k)] 



(BIO) 



where is the Fermi-Dirac distribution function and /9^^(/co,k) is the free spectral function for the fermion field. 



px, (fco, k) = 27r (- + mxj + (k)) 



(Bll) 



with ojXi (k) = ^/k' 



As in the previous decay mode case, when (jBlip is used in (jBlOp , we can again identify decay and Landau damping 
terms. The momentum integration is carried out also in a very similar way, giving as the final result, from the 
definition Eq. (j6.2p , that the decay width of a scalar boson X decaying into a fermion-antifermion pair Xi + X2 is 
given by 
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87riLJx(p) 
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/ \l + e 3 
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l + e- 
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l + e 
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{-pI + P^ + {mxi - mx2 



(B12) 



3. Decay of a fermion to a scalar boson and a fermion 

The last case is the decay of a fermion X to another fermion Xi and a scalar boson X2- The Yukawa coupling 
among the fields is again gx ■ The imaginary part of the fermion field X self-energy is in this case given by 

/d^k 
j^-^npipo - ko)pxi{po - ko,p -k)nB{ko)px2{ko,'^) ■ (B13) 

Using the expressions for the free spectral functions, Eqs. (jB2l) and (jBllI) . considering po > 0, the branch cuts are 
still given by -\m\^ - m\^\ < pi - < {mx^ - mxa)^ and (mxi + "i-Xa)^ < pI ~ < oo- We obtain 
ImSix = ImS^'^^'' + ImS^^ , where the decay term is 

ImE^^'=''''(po,p) = 27r5| [ y^-^l 1 T^i ttt {7"bo - (k)] - 7-P + 7-k + } 

^'^ J (27r)'^ 4wxi(p- k)tJx2(k) ^ 

X [l + ni3(wx2(k)) -nF(po-wxi(k))](5(po-a;xi(p-k) -wxjk)) , (B14) 
and the Landau damping term is 

/cPk 1 
TT^'a 1 r\ 7r^{-^°[Po-^^xA^)] +7-P-7-k-mxi} 
{2ttY Aujxi (P - k)a;x2 (k) 

X [nF(po + '^Jfi(k))-l-nB(a;x2(k))](5(po-wxi(p-k)+wjf2(k)) . (B15) 
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The momentum integrals in the above expressions again foUow analogously to the scalar boson case. This fermion 
self-energy can now be expressed in the form Eq. (|6.4p . from which explicit expressions for ckq, cxp, and 5^ can be 
identified. These will be listed below. 

The decay width in this case can be determined from the poles of the dressed fermion propagator, + mu x + 
«ImSx)~^, given by 



± ujx{p) ± i 



LOx{v) i&p + ao) H 7-T {am - ap) 



wx(p) 

where the decay width of fermion X to fermion Xi and boson X2 is found to be 



±ujx{p) ± iTx{po,p) , 



(B16) 
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wx(p) 

The explicit expressions for the a coefficients are 
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p2 + p2 _ rn2^ + r„2^ „o 

1 _ e— +/^ 



p^ - (mxi + mx^)"^] 



S^IpI I 2p2 

(P§ + P^ - mji, + - 2poW-) 



In 



2p2 



In 



1 + g-\Po-^+\/T 



1 -|_ g-bo-"-l/T 



9[pi - p' - {mx, +7nx^Y] 



9'xT'po 
87r|p|3 

Li2 

9%T 
8^|p| 



{Li2 (e-"+/^)-Li2 (e--/^)+Li2 
} 9[pI - P^ - + mx^f 



-\po-u+\/T 



3^ + p2 - m^^^ + mj^^ - 2poUj+) 



2p2 



In 



1 - 



/T 



+ p2 - m^^ + TO^^ - 2poU}-) 



2p2 



In 



1 + p-(po+w+)/T 



1 + e-(Po+"-)/'^ 



-Po + P +{mxr -mx^y] 



^ {Li, -Li, (.-"-'-J + Li, 

} ^[--Po + P^ + {mx^ - mx^f] , 



_g-(po+i^+)/T 



-(po+tj_)/T 



Q!m = „ , , u)-)0[pQ--p -{mx^+mx^}] 



SttIpI niR^x 
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9x ^j.^ , 1^ 



9x ^^J^ , 1^ 



1 - e 



-UJ+/T 



SttIpI niR^x 



1 + e~IPo~"+l/'^ 
1 - e-'^+Z^ 



1 + e-(po+<^+)/'3 



where Li2(x) is the dilogarithm function, 



and ijj± are given again by Eq. 



hi 
In 



1 - e 



-u]-/T 



1 + e-IPo-'-'-l/T' 
1 + e-(po+"-)/T 



6'[Po -P^ - {mx,+mx^f] 



(B20) 



(B21) 
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